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SOME POSSIBILITIES OF USING GAS MIXTURES OTHER THAN AIR 

IN AERODYNAMIC RESEARCH 1 

By_lJF.A.\\R.T!HA PMA X 


SUMMARY 

a l study is made of the advantages that can be realized in 
compressible-flow research by employing a substitute heavy gas 
in place of air. A lost heavy gases considered in previous in- 
vestigations are either toxic, chemically active, or (as in the case 
of the Freons) have a ratio of specific heats greatly different 
from a ir. The present report is based on the idea that by properly 
mixing a heavy monatomic gas with a suitable heavy polyatomic 
gas, it is possible to obtain a heavy gas mixture which has the 
correct ratio of specific heats and which is nontoxic, nonflamma- 
ble, thermally stable, chemically inert, and comprised of com- 
mercially available components. 

Calculations were made of wind-tunnel characteristics for 
(>3 gas pairs comprising 21 different, polyatomic gases properly 
mixed with each of three monatomic gases (argon, krypton, and 
xenon). For a given, Mach number, Reynolds number, and 
tunnel pressure, a gas-mixture wind tunnel having the same 
specific-heat ratio as air would be appreciably smaller and would 
require much less power than a corresponding air wind tunnel, 
x inalogous though different advantages can be realized in com- 
pressor research and in firing-range research. 

'J'hc most significant applications, perhaps, arise through 
selecting and proportioning a gas mixture so as to have at 
ordinary wind-tunnel temperatures certain dimensionless char- 
acteristics which air at flight temperatures possesses but which 
air at ordinary wind-tunnel temperatures does not possess. 
Characteristics which involve the relaxation time (or bulk vis- 
cosity), the variation of viscosity with temperature, and the 
variation of specific heat with temperature fall within this 
category. Other applications arise in heat-transfer research 
since certain gas mixtures can be concocted to have any Prandtl 
number in the range at least between 0.2 and 0.8. 

INTRODUCTION 

'J'lic reasons for considering gases other than air as possible 
lest media for compressible-flow research stem primarily 
from the. relatively low speed of sound in certain gases. In 
general, the heavier the gas, the lower the speed of sound at 
a fixed temperature. Hence, experiments conducted at a 
given Mach number in a heavy gas will be conducted at 
lower velocity than the corresponding experiment in air. 
Some significant advantages of conducting wind-tunnel, 
firing-range, and compressor experiments at a reduced 
velocity have been noted previously in the work of Tlieo- 


dorsen and Regier (ref. 1). Smelt (ref. 2), Kantrowitz (ref. 
3), Huber (ref. 4), Buell (ref. 5), Donaldson (ref. 6). and 
von DoenhofF and Braslow (ref. 7). Some uses of gases 
other than air in shock tube research are discussed by Duff 
in reference 8. These investigations show that the substi- 
tution of a heavy gas for air offers the possibility of: (1) 
extending the range of existing research apparatus; (2) 
achieving greater economy of construction and operation of 
large high-speed wind tunnels; and (3) providing greater 
facility in obtaining data for special types of research. 
These advantages will be discussed briefly in the order listed. 

The possibility of extending the range of existing apparatus 
was clearly demonstrated by the experiments of Theodorsen 
and Regier, which appear to be the first compressible-flow 
experiments in which supersonic aerodynamic data were 
obtained in a gas other than air. By rotating propellers in 
Freon 113 (CCUFCClFo, having a speed of sound 0.39 times 
that of air, and a density 6.5 times as groat) the}- were able 
to achieve tip Mach numbers of 2.7; whereas the highest 
tip Mach number achieved with the same apparatus using 
air was 1.0. Also, much higher Reynolds numbers were 
obtained. In a similar fashion, Buell employed the Freons, 
and Donaldson employed xenon as the test medium through 
which projectiles were fired in order to obtain Mach numbers 
much higher than could be obtained with the same apparatus 
by firing through air. An extended range of operation of a 
wind tunnel has been demonstrated by von Doenhoff and 
Braslow who report that the maximum attainable test- 
section Mach number of the. Langley low-turbulence pressure 
tunnel was increased from 0.4 to 1.2 by replacing air with 
Freon 12. 

A possibility of achieving greater economy by employing 
a substitute heavy gas in a wind tunnel arises because the 
test-section dimensions are smaller and the horsepower 
required is only a fraction of that required for an air wind 
tunnel operating at the same Mach number, Reynolds 
number, and pressure. Smelt has made a general study 
pointing out this possibility for various inert gases that have, 
unfortunately, a specific-heat ratio different from air. The 
importance of reduced power requirements, if accomplished 
without sacrifice in ratio of specific heats, needs no elabora- 
tion in view of the power of modern wind tunnels. 

The third advantage mentioned, namely, facilitating 
certain types of aerodynamic research, arises in two separate 
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ways. First, for a given Macli and Reynolds number, 
the use of a heavy gas with its low velocity of sound enables 
experiments to be conducted at lower velocity and lower 
pressure than if air were employed; this results, for example, 
in lower stresses in wind-tunnel models as well as in models 
launched from a gun. Similarly, the reduced velocity 
enables experiments on the release, of objects from super- 
sonic vehicles to be conducted at smaller scale while simu- 
lating the required value of the Fronde number (£'» 2 /gl)- 
Second, and more important, the use of certain gas mixtures 
other than ah - enables some of the dimensionless parameters 
pertinent to hypersonic flight through air to be simulated at 
the usual wind-tunnel temperatures whieli cannot be simu- 
lated by use of air at these wind-tunnel temperatures. 
Several of such possibilities are discussed later in this report. 

As regards the various disadvantages involved by using 
substitute hea\ T 3 r gases in aerodynamic research, the main 
ones are: ( 1 ) The ratio of specific heats ( 7 ) differs from 
that of ah- for the heavy gases proposed thus far; ( 2 ) the 
use of any gas other than air results in more severe practical 
problems of operating a research facility'.- A general 
appraisement of the various advantages relative to dis- 
advantages cannot be made since this would depend on each 
particular case. The purposes of this report are to study 
various heavy gas mixtures which appear satisfactoiy for 
aerodynamic testing and to evaluate the magnitude of the 
advantages such gases offer. Knowledge of these advantages 
will enable a comparison with the disadvantages to be made 
in a given case. 

. The present research was begun upon conception of the 
following simple idea : 3 The proper ratio of specific heats 
( 7 = 1 . 4 ) can be achieved with a heavy gas by suitably 
mixing a heavy monatomic gas (yI>lA) with a heavy 
polyatomic gas (y<(1.4). Exploitation of this simple idea 
turned out to be unexpectedly arduous because of difficulties 
involved in obtaining adequate chemical, physical, pharma- 
cological, and thermodynamic data on the many polyatomic 
gases known to modern chemistry. Primarily as a result 
of recent developments in the field of fluorochemicnls, it was 
possible to find among the known polyatomic gases over :30 
which were indicated by available data to be sufficiently 
nontoxic, nonflammable, chemically inert, and thermally 
stable to justify consideration. Evaluation of the relative 
advantages of each polyatomic gas for aerodynamic testing 
requires calculating the thermodynamic properties from 
spectroscopic data, determining the proper proportions for 
mixing from thermodynamic properties, and then computing 
the density and viscosity that would result, when mixed 
with each of three monatomic gases considered (argon, 
krypton, and xenon). 

2 A possible disadvantage sometimes expressed (c. g., ref. 9) is that the relaxation time for 
molecular vibrations in certain polyatomic gases would invalidate data obtained at super- 
sonic velocities. However, for most gases considered here, this is believed not to be the case. 
Reasons for this are discussed later. 

3 After completion of the main body of calculations, short unpublished note by E. F. Relf, 
in England, became available Jn which the possibility of obtaining the proper ratio of specific 
heats through mixing gases also was pointed out. Itelf made a rough estimation for a mixture 
of sulfur hexafluoride and xenon, which is one of the mixtures considered in flic present 
research. 


NOTATION 

a speed of sound 
c speed of light 
<*,, specific heat per unit mass 
C p specific heat per mole 
('ll wind-tunnel compression ratio 
<j acceleration of gravity 
li Planck’s constant 
H molar enthalpy 
TIP horsepower 
A- Boltzman’s constant 
l characteristic length of model 

L characteristic dimension of wind-tunnel lest section 
m molecular weight (29 grams per mole for air) 

.1/ Mach number 
v number of atoms per molecule 
p pressure 

, 1 . 

<[ dynamic pressure, -^pu- 
ll universal gas constant per mole (1.987 cal. mole -1 
0 K -1 ) 

Re Reynolds number 
T temperature 
u velocity 
V volume of gas 

w mass-flow rate, (pi/) times (cross-section area) 
x mole fraction, equal to ratio of partial pressure to 
mixture static pressure and equal to fraction by 
volume 

t relaxation time 
7 ratio of specific heats 
k bulk viscosity coefficient 
p mass density 
p viscosity coefficient 

w wave number (cm -1 , from spectroscopic data) 

SUBSCRIPTS 

00 free-steam conditions in test section 

c critical conditions for gaseous phase of a compound 
t total conditions for gas brought isen tropically to rest 
(wind-tunnel reservoir conditions) 

° 0° centigrade 

1 monatomic gas 

2 polyatomic gas 

SUPERSCRIPTS 

- quantity divided by corresponding quantity for air 

METHODS OF SELECTING GASES AND COMPUTING CHAR- 
ACTERISTICS OF VARIOUS GAS MIXTURES 

SELECTION' OF GASES 

The selection of monatomic gases for consideration is 
simple, since only four such gases are known — argon, 
krypton, xenon, and radon — that have a molecular weight 
greater than air. Radon, the heavist of all, unfortunately 
must be excluded because of its radioactivity. The remaining 
three gases arc completely inert. Some of their properties 
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pertinent lo the present analj’sis are listed in table I. Al- 
though argon is the only one now commercially available in 
large quantities, both krypton and xenon will be considered 
since (hey might be of use for special research in small 
apparatus, and since their availability undoubtedly will 
increase in the future. 4 

'Pile selection of polyatomic, gases for consideration was 
rather tedious. Of the approximately 13,000 polyatomic 
compounds listed in the 1952 Handbook of Chemistry and 
Physics (ref. 12), 103 were found to be gases at 0° C. Com- 
pounds that arc not gaseous at 0° C were not considered. 
Most of the gases were not considered further because they 
possessed some intolerable characteristic such as being 
highly poisonous (e. g., CO Cl-., >S 2 F 10 , C 2 X 2 , S0 2 ), toxic (e. g., 
(TT 3 F, CF 2 Br 2 , C'H 3 CH 2 F), combustible (e. g., COS, all 
hydrocarbon gases), unstable (e. g., CHCC1, C1X 3 , PSF 3 ), 
or generally active chemically (e. g., BF 3 , FT), CF 3 SF 3 , 
CF 3 OF). Still others were eliminated for less obvious 
reasons, such as reacting with water vapor to yield corrosive 
products (e. g., SiF 4 , PF 3 , COF 2 , TcF 6 ), or having a tendency 
to polymerize under pressure (e. g., CC1FCF L , CH.CF-, 
CF 2 CF 2 ), or being active when in the presence of mercury 
(F 3 CXNCF 3 ), or imparting, upon inhalation, an offensive 
odor to the breath which persists long after inhalation 
(TeF 0 ). This elimination process left 16 polyatomic gases 
for further consideration. It was noticed that all but two 
((XX and NIO, the lightest two on the list) contained at least- 
two atoms of fluorine per molecule. Because fluorine 
appeared Lo be the key to the problem, recent technical 
publications on fluorine chemistry were scanned (especially 
refs. 13, 14, and 15) for further information. From this 
search 15 additional polyatomic gases were found which 
met the required conditions of being nontoxic, noncom- 
bustiblc, chemically inactive, and thermally stable. Among 
these 15 recently synthesized gases were several (e. g., 
C'BrF 3 ) which turned out to be the most efficient for wind- 
tunnel use. Undoubtedly, additional new polyatomic gases 
will continue to be synthesized in the future inasmuch as a 
large industrial effort is being expended on fluorine chemistry. 
Quite a few apparently satisfactory gases (c. g., C 3 C1F 7 and 
(311 Fr.) are known to be omitted from the list. Conse- 
quently, the present report which primarily indicates current 
possibilities of gas mixtures in reducing wind-tunnel power 
is by no means exhaustive and might considerably under- 
estimate possibilities existing in the future. 

It is pertinent to focus attention here on the extreme prop- 
erties of fluorochemicals. Their chemistry often is totally 
unlike that of analogous halide compounds. Fluorochemicals 
include some of the most toxic compounds known, yet also 
include the most inactive and nontoxic compounds known. 
Them thermal stability and resistance to electric breakdown 
are unusually high. Most of these extreme properties can 
be illustrated by a few examples. First consider the follow- 
ing two gases: chloroform, CHC1 3 , and fluoroform, CHF 3 . 

* references 10 and 1 1 for information on the possible future availability of the monatomic 

Rasas. 


The toxicity of chloroform is well known. Fluoroform, how- 
ever, is indicated to be completely inert physiologically (ref. 
16). Guinea pigs can live in an atmosphere of 80 percent 
fluoroform and 20 percent oxygen without being affected 
either during the test or afterward. Sulfur hexafluoride. SF 6 , 
also is indicated to be physiologically inert on the basis of 
similar tests with small animals (ref. 17). On the other hand, 
the deceptively similar sulfur pentafluoride, (SF 5 )_>, is ex- 
tremely poisonous, since, for concentrations as low as one 
part per million — -in which concentration phosgene is harm- 
less — -sulfur pentafluoride is lethal (ref, IS). Further illus- 
tration of the extreme properties of fluorochemicals is 
provided by carbon tetrafluoride, CF 4 . The thermal decom- 
position temperature of CF 4 has not yet accurately been 
determined since it is so high as to be in the vicinity of the 
carbon arc temperature (ref. 19). Carbon tetrafluoride has 
been heated in the presence of many metals up to the 
temperature at which glass softens without any reaction 
occurring. (See ref. 14, p. 433.) Sulfur hexafluoride is 
equally remarkable in its chemical inactivity at all tempera- 
tures below the softening point of glass. 

The chemical symbol, molecular weight, chemical name, 
and trade name of 31 polyatomic gases considered as possible 
components for aerodynamic test media are listed in the first 
four columns of table II. A symbol designating the general 
classification of each gas is listed in the fifth column. The 
meaning of each symbol is explained in the footnote of the 
table. Listed in table III are certain data, to be utilized 
later, regarding the boiling point, critical constants, viscosity 
coefficient, and spectroscopic frequencies of the various poly- 
atomic gases. The sources of data, for each gas are indicated 
by reference numbers in the last three columns of table II. 
These data were compiled from references 13 to 16 and 18 
to 56. 

computation of thermodynamic properties 

The first step in evaluating aerodynamic qualities of any 
gas is to determine its thermodynamic properties such as 
specific heat and enthalpy. When two gases are mixed, the 
resulting thermodynamic properties can be calculated from 
the thermodynamic properties of the individual components. 
Using subscript 1 to denote the monatomic component, and 
subscript 2 lo denote the polyatomic component, the follow- 
ing equations apply for binary mixtures of thermally perfect 
gases: 

molecular weight m=riin l -r :r 2 rrt 2 (1) 

equation of state p=pRTjm (2) 

specific heat C v -]rR=C p —XiC Pi -rx 2 C p „ (3) 


X\C p . -pXoUp, 


ratio of spocific l.oate y~^~Z^)+^=R) 
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where R is the universal gas constant (1.9S7 eal. °Iv _I 
mole -1 ), and x denotes the mole fraction of a component 
(equal to the fraction by volume, or the fractional partial 
pressure of the component) and is subject to the relation 
2 i+a- 2 = 1. For all monatomic gases, 

/^4.97 cal. o K -i mo ] 0 -i 


Hence, from equation (4), it follows that- for 7=7/5, 


1 — a*i=ar ; 


R 2 ) 


(6) 


The thermodynamic properties of a polyatomic gas can he 
accurately computed from infrared and/or Raman spectro- 
scopic data, using equations developed by methods of statis- 
tical mechanics (see. ref. 57. for example). For linear 
molecules (e. g., C0 2 and X 2 0), wherein the atoms are ar- 
ranged in a straight line, 



The various wave numbers w,- determined from spectroscopic 
data for each gas are listed in table III. The constant 
hc/k is equal to 1.4385 °Iv cm. For nonlinear molecules not 
subject to internal rotation (e. g., CF 4 , CHF 3 , SF s , etc.) 



For nonlinear molecules subject to a torsional mode of oscilla- 
tion (e. g., CF 3 CF 3) CHFoCCIFo, and all other ethane-like 
molecules) the above summations were employed only for 
the 3n-5 normal modes, and the torsional contribution was 
computed from the Pitzer-Gwinn tables (reproduced in 
ref. 5S) using values for the potential barrier estimated from 
spectroscopic data. Computations were greatly facilitated 
through use of the tables of reference 58. Typical curves 
showing the molar specific heat as a function of temperature 
for several of the various gases considered are presented in 
figure 1. Curves of this type, together with equation (4), 
were used for each gas to determine the. mixture proportions 
necessary to obtain a value of 1.4 for 7. Similar data for 
the enthalpy were used to compute Mach number, pressure 
ratio, etc., for the different flow conditions considered. 


It is noted that the thermodynamic properties calculated 
from spectroscopic data agree well with direct calorimetric 
measurements. This is illustrated in the following table: 


Gas 

T, Q C 

Calculated C r 
(from spectro- 
scopic data of 
table III) 

Experimental 
C v (from direct 
calorimetric 
measurements) 

co> 


7.71 

1 

7.78] 


20 

, 8. $2 

S. Sol-ref. 57 


220 

10.62 

10. 63) 

CC1/I-Y 

—30 

15.6 

15.7 ] 


0 

1G.G 

1G.7 

18.0 | ref - 01 


45 

17.8 


90 

IS. 9 

19.1 1 

CI-'.CI-Y 

-50 

21.1 

21.0 j 


10 

24.7 

21. -1 |rcl. 22 


90 

i 2S.7 

j 

2S.2 ] 


COMPUTATIONS OF VISCOSITY 

It will be seen subsequently that for a given Reynolds 
number and Mach number, the power required by a wind 
tunnel is directly dependent on the gas viscosity coefficient. 
The viscosity of a mixture of two gases can be approximately 
calculated from the equation (see ref. 59): 



It is noted that a simple linear interpolation between the 
values of ju for each component is not an adequate method 
of estimating the viscosit}" of the mixture. This is illustrated 
by several examples in figure 2. Curves of this type were 
constructed for all gas mixtures considered. 

Experimental data are available for the viscosity mi of the 
three monatomic gases, but very little data arc available for 
the various polyatomic gases. Fortunately, the coefficient 
of viscosity can be estimated quite satisfactorily from 
knowledge, of the critical temperature T c , critical volume 
V c , and molecular weight- by the equation 



This equation may be deduced from the general equations 
for gases having intcrmolccular forces dependent only on 
powers of the distance between molecules. (Sec. refs. 59 
and 60). Rather than to use the elaborate theoretical 
evaluation of the function F{T/T C ), the more simple empirical 
evaluation of reference 61 is used: 


F 


27.1 TjT c for 

©H for 7’> 1 .6 T r 

1+0.0952 (jr-l.QJ 


r 02a.) 


where T is in °K, m is in gm./mole, V c is in ec/gm. mole, and 
n is in micropoise. A comparison of the calculated and 
experimental viscosity coefficients of all gases considered for 
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which experimental data could be found is given in the 
following table: 



Viscosity at 0° C 


in rmcropoises 


Culcu- 

Expert- 1 


| hited 

mental , 

Ar 

1 210 

1 210 

Kr 

' 229 

! 230 

Ne 

211 

210 

SF 8 

146 

142 

CCIjFj 

. 114 

118 

CHCIFj 

1 119 

120 

CO,.... 

134 

137 

NjO 

132 

135 

Air.., 

1 170 

172 


This agreement is satisfactory. Consequently, in the calcu- 
lations of wind-tunnel power requirements, experimental 
values of n are employed for the 9 gases listed above, but 
values computed from equation (12) are employed for the 
remaining gases for which experimental data are not avail- 
able. .Since molecules of the above 9 gases are either simple 


or approximately spherical in structure, the accuracy of 
calculations from equation (12) is expected to be less for 
greatly different molecular structures, such as that of 
CjFio, but probably is sufficient for present purposes. 

WIND-TUNNEL POWER AND SIZE 

The general equation for compressor power required to 
maintain a steady rate of mass flow a; of a perfect gas through 
a compression ratio CR starting with an initial temperature 
T, is given by the equation 

Hp=J V'C I T i ^CEyy-i)n-i\ (13) 

v 

where J is a numerical constant depending only on the units 
employed, v is the combined adiabatic and mechanical 
efficiency (assumed to be independent of the gas), and c p is 
the specific beat per unit mass. It is noted that c p in equa- 
tion (13) is related to the molar specific heat C p and the 
molar universal gas constant R through the equation 

mc p =C p =~^ R 


(14) 
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Figuke 2. — Some typical examples of viscosity coefficient for gas 

mixtures. 


By using a bar to denote a quantity divided by the corre- 
sponding quantity for air (e. g., m = m/(m) ; , ir , w=w/{w) air , 
HP =HP j{HP) ai „ etc.), there results 

HP== T,Y (y, T[Y(y, CP) (15) 

m in 

where Y(y,CR)=~-^ [(C7£) (Y-I)/7 — 1], and L is a character- 
istic dimension of the test section. A more convenient 
equation is obtained by introducing the Reynolds number 

T t Y( 7, CR) (1 6) 

m 


from which it is apparent that, as Smelt (ref. 2) lias pointed 
out, power economy can be achieved by employing a gas 
with low viscosity and high molecular weight, or by testing 
at the lowest temperature possible, or by testing in the small- 
est facility (highest operating pressure) that will yield the 
given Reynolds number. The quantities Re and L are rela ted 
to the relative operating tunnel pressure p ( 




or, finally, 


where 


7 _ 


/(T-I/Js 



( 17 ) 


Equation (17) shows that by increasing the molecular weight, 
or decreasing the viscosity, the size of wind tunnel required to 
obtain a given Reynolds number can be reduced. Hence, 
an alternate equation for relative wind-tunnel power is 


Y (y,CR) 

P'&T* 


(is) 


Equations (16), (17), and (18) are a generalization of similar 
equations developed by Smelt. 

For subsequent computations, the above equations can be 
simplified. It will be assumed that CR depends only on M„. 
(To achieve this in an air tunnel, the compressor speed would 
have to be reduced for a heavy gas.) For the case of tran- 
sonic wind tunnels CR is the order of 1.3_or less, and direct 
numerical calculations show that_0.97<^F(7,1.3)-O for any 
7 between 1.1 and 1.4. Hence, T-=l is a good approxima- 
tion for transonic wind tunnels irrespective of 7. For the 
case of supersonic wind tunnels, only_ gas mixtures having 
7 — 1.4 arc considered, for which TQ_-4,C7?) — 1. Conse- 
quently, in all cases of concern here, T=1 is an adequate, 
approximation. A similar argument shows that the ratio 
1^(7, Cf?)/_/ (7, il/oo) is approximately unity for the range of 7 
and considered here. _Inasmuch as all calculat ions will 
be made on the basis that T, — 1, the above equations simplify 
to 

(16a) 

m 


R n 
Prjm 


(ipfiny 

IhCmf 1 


(17a) 


11 8a) 


Four illustrative cases will be considered: (1) given Re, 
and L; (2) given d/„, L, and p,; (3) given M„, L, and HP; 
and (4) given d / Re, and p t . Since for all cases T, = 1 , 
the reduction in velocity is always u m = (m)~ ii . 

Case (1) : Given Re, and L , — This case corresponds 
to operating a given wind tunnel (L — 1) with a gas other than 
air at the same values of d/„ and Re as for air. The horse- 
power required for the gas mixture relative to that for air is 
HP—Jipn. However, from equation (17a), it is seen that 
M for given values of Re and L. Hence, hi order 
to achieve the same Re with a heavy gas as with air, the 
given tunnel would be operated at a lower pressure, resulting 
in model loads that arc proportionately lower, since q—pt= 
m/VTF 

Case (2): Given M L, and p t . — This case corresponds 
to substituting a gas other than air hi a given whid tunnel and 
operating at the same pressure and Mach number as for air. 
From equation (17a) it follows that This indi- 
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oat os that, higher Reynolds numbers would be achieved with 
a heavy gas than with air, while (according to eq. (16a)) 
simultaneously drawing less power by the ratio II P — (m) ~ ’ rl 
Since p t — 1, model stresses are unchanged for this case. 

Case (3): Given L, and HP. — This case corresponds 
to replacing air in a given tunnel by a gas mixture and draw- 
ing the same power for a fixed Mach number. The Reynolds 
number, from equation (16a), would be increased by the 
factor In order to draw the same power as with 

air, the pressure and, hence, model loads would have to be 
increased by the factor p, = y7w. 

Case (4): Given Re, and p t . — This case corresponds 
to conditions that might be prescribed in the initial design 
of a new wind tunnel whose size is to be determined by the 
required values of A/„ and He, by the selected value of p t , 
and by the gas employed. From equation (18a) it is evident 
that in this ease, the horsepower required for a heavy gas 
mixture is less than for air by the factor 

(l9) 

This represents a greater power reduction than in the cases 
above because the gas-mixture wind tunnel, according to 
equation (17b), would be smaller than the corresponding air 
tunnel by the factor 

I=HAm (20) 

The gas-mixture wind tunnel would produce the same aero- 
dynamic data as the larger air tunnel with the same model 
stresses. Computations based on these latter two equations 
are presented later. 

MACH NUMBER LIMIT FOR CONDENSATION 

Although any of the polyatomic gases would be usable in 
apparatus where low temperatures are not involved, such 
as in a subsonic wind tunnel, a firing range, or a compressor 
research apparatus, only a limited number would be useful 



Figure 3. — Method of determining saturation limit. 
377737— 3<i— 2 


in a supersonic wind tunnel where the static temperatures 
encountered are low and can result in condensation. An 
essential step in evaluating various polyatomic gases, there- 
fore, is to determine the approximate useful Mach number 
range for each gas. The method of determining this was 
to draw the isentropic expansion curve for the partial pres- 
sure of each polyatomic gas on log-log paper together with 
the curve for saturation vapor pressure of that gas. As 
indicated in figure 3, the intersection point yields T min , from 
which the maximum Mach number can be calculated from 
the relation 


y \RT RTJ 


where 7 and II/T are determined from spectroscopic 
data by equations (4) through (10). Graphical solutions of 
this type were made for operating conditions of approximately 
40° C (100° F) total temperature and several atmospheres 
total pressure. These solutions enabled the appropriate 
polyatomic gases to be selected for each of several design 
Mach numbers considered in the evaluation. It is emphasized 
that the graphical method included consideration of the vari- 
ation of specific heat with temperature. As will be seen sub- 
sequently, such variation is of dominant importance in 
determining the limiting Mach number to which certain 
gases can be expanded without reaching saturation. Simpli- 
fied criteria, such as boiling temperature, are quite inadequate 
for determining the saturation limit. 


RESULTS AND DISCUSSION 

WIND-TUNNEL POWER AND SIZE 

Complete computations could not be made for 10 of the 
31 polyatomic gases listed in table II, since the necessary 
spectroscopic data and critical constants were not available 
for these 10 gases. Calculations have been made, however, 
of the relative wind-tunnel power requirements for 63 gas 
pairs comprising mixtures of the 21 polyatomic gases for 
which data are available with each of the 3 monatomic 
gases. Mixture proportions were determined from equation 
(6) by substituting a value of 1.4 for y a . The calculations 
were made for three design Mach numbers, 1.3, 2.5, and 3.5. 
Results are tabulated in tables IV (a), IV (b), and IV (c), 
respectively. It should be noted that all 21 polyatomic gases 
are included in table IV (a), since all are free from condensa- 
tion at M„~ 1.3. Only 7 gases for M a — 2.5 (table IV (b)) 
and 2 gases for AG, = 3.5 (table IV (c)) have sufficiently low 
boiling points to avoid saturation for t ho assumed condition 
of 40° C total temperature. Operation at Mach numbers 
above about 3.5 would require total temperatures higher than 
the value of 40° C arbitrarily assumed. 

The values of relative horsepower HP given in table IV 
represent case (4) mentioned previously. This case corre- 
sponds to a comparison with an air wind tunnel for the same 
Af„, Re, and p t . If desired, any of the other three cases 
mentioned, or any of the individual relative quantities such 
as or ~p a , can be readily calculated from the values of 
7 n and Ji listed in the tables and from the appropriate equa- 
tions developed earlier. 

The tabulated results show that there is no single poly- 
atomic gas which is best for all Mach numbers or for use with 
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all monatomic gases. For each monatomic gas, though, there 
are several polyatomic gases which appear roughly equal in 
their ability to require low wind-tunnel power. Of the gases 
that currently are commercially available and not classified 
questionable in table II, the most efficient ones for use with 
argon are: CH 2 F 2 , CBrF 3 , and CCL>F 2 for M m — 1.3; CBrFa 
and SF 6 for M a = 2.5; and CF 4 for A/*, =3.5. In general, 
the horsepower and size of gas-mixture wind tunnels relative 
to equivalent air wind tunnels are very approximately as 
follows: 


1 

" 

Arpon 

Mixtures 

Krypton 

::j 

Xenon I 

1 

Power relative to air. U1 J 

0.3 

0.2 

0. 1 

Test-section dimension relative loair. 




L 

0,7 

0.0 

0 5 1 


T hese figures indicate significant advantages of gas mixtures 
as wind-tunnel test media. 

It is interesting to note that the heaviest polyatomic 
gas — C 4 Fio, having 7%— 8.2 — does not result in either the 
heaviest mixture or the best mixture This is due to the 
high molar specific heat of CiF 10 which requires that only a 
small amount (6.5 percent) be mixed with a monatomic gas 
in order to obtain 7 = 1 . 4 . However, for applications to 
low-speed subsonic research where the value of 7 does not 
matter and, hence, a mixture is not needed, the heaviest 
gas is the best. For example, pure C 4 Fi 0 , having 7=1.06, 
would require only 1.8 percent of the power, and an appa- 
ratus one-fourth the size of that required by air for a given 
, Be, and p, In this respect pure C 4 F 10 is several times 
more efficient than pure Freon 12. (See table III where 
corresponding values are given for all pure polyatomic gases 
considered.) It is evident that pure polyatomic gases hav- 
ing 7 ~ 1 .1 are more efficient in reducing wind-tunnel power 
requirements than are gas mixtures having 7 = 1 . 4 . 

A possibility that should not be overlooked is that a 
value of 7 close to 1.4 may not be necessary for transonic 
wind tunnels or other transonic research apparatus. I 11 
view of the surprisingly small differences observed at tran- 
sonic speeds by von Doenhoff and Braslow (ref. 7 ) between 
uncorrected results in Freon 12 (7=1.13) and in air ( 7 = 1 . 4 ), 
it would appear that a reduced value of 7 , perhaps between 
1.3 and 1.2, might yield transonic data directly applicable 
to air for practical purposes. If a considerably reduced 
value of 7 is satisfactory in transonic wind tunnels (such is 
not anticipated in supersonic wind tunnels), then the possi- 
ble advantages of gas mixtures for transonic speeds are 
greater than the above calculations for 7 = 1.4 would indi- 
cate. The following table for CBrF 3 -Ar mixtures at 
A/ co = 1.3 illustrates this: 


Case (4): given 
Re, pt. and A/co 


y 

IIP 

L 

1 .4 

0 37 

0. 7-1 

1.35 

.30 

. 00 

1.3 

.22 

.02 

1. 25 

. 10 

.53 

I. 2 

. 10 

. 44 

1. 15 

-00 

.37 


Case (3): 
given 
UP, A/co, 
and L 

~Rc 

2.0 (from (able IV(a)i 

2.3 
2.8 

3.4 

4. 5 

0.2 (pure CBrFa 


These figures, which are typical of many of the argon mix- 
tures, show that by reducing 7 to about 1 . 28 , the power 
requirement for an argon gas mixture at transonic speeds 
would be about one-half that for the same type mixture 
proportioned to yield a 7 of 1 . 4 . This corresponds to one- 
fifth the power of an equivalent air wind tunnel operating 
at the same Reynolds number, total pressure, and Mach 
number (case (4)). Also, the Reynolds number of an exist- 
ing transonic wind tunnel having a fixed horsepower (case (3)) 
would be about tripled if air were replaced by an argon gas 
mixture having 7= 1 . 28 , but only doubled by a mixture hav- 
ing 7=1.4. The advantages of reducing 7 to the lowest 
practical value are apparent from the table. For xenon 
mixtures, however, there is less to gain by reducing 7 below 
1 . 4 . It would be desirable to conduct some experiments 
varying the proportions of a gas mixture in order to deter- 
mine to what limit 7 can be reduced in a transonic flow and 
still yield data directly applicable to air for practical purposes. 

DEVIATIONS from thermal and caloric perfection 

Under wind-tunnel conditions air generally is considered 
as being both thermally perfect ( pm — pBT ) and caloricallv 
perfect ( 7 =constant). Most other gases can not be so con- 
sidered. Calculations of the degree of deviation from 
thermal perfection for the various argon gas mixtures are 
presented in Appendix A. It will suffice here to state that 
such deviations are indicated to be unimportant for normal 
operating conditions of wind tunnels. The smallness of 
these deviations is due in part to the monatomic (argon) 
component, which by itself is almost thermally perfect, 
and, in part, to the fact that the polyatomic component is a 
fluorochemical. The unusually low intormolccular forces in 
fluorocliemical gases result in small deviations from thermal 
perfection. 

Deviations from caloric perfection, on the other hand, are 
not small and represent an important technical consideration. 
The specific heat of polyatomic gases varies widely, as figure 
1 clearly indicates. Thus, a gas mixture proportioned to 
yield 7 = 1.4 for one Mach number will not yield the same 
value of 7 at other Mach numbers. This may limit the 
Mach number range over which a single gas mixture could 
be used satisfactorily. Also, at a fixed Mach number, 7 
will vary with temperature. Under wind-tunnel conditions, 
7 for air is quite constant, but varies under flight conditions 
where higher temperatures are encountered. 

To aid in evaluating the importance cf caloric imperfec- 
tions, the variation of 7 with temperature has been computed 
for various gas mixtures and for air under both wind-tunnel 
and flight conditions. In these computations the spectro- 
scopic frequencies listed in table III wore employed to 
determine 7 for the gas mixtures, whereas the tables of 
reference 62 were employed for air. Instead of showing the 
computed variation as a function of some temperature 
parameter, an enthalpy parameter is em- 

ployed which always is zero for free-stream conditions and 
always is unity for reservoir conditions, regardless of the 
nature of the gas. Some of the results are shown in figure. 
4. Air under flight conditions at a given M«, is represented 
by a cross-sectioned band because ambient temperatures in 
the atmosphere vary considerably with altitude. Gases 
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(a) Gas mixture of 0.25 CBrFj with 0.75 monatomic. 

(b) Gas mixture of 0.4 CFi with 0.6 monatomic. 


Figure 4. — Variation of specific-heat ratio between free-stream and 
stagnation conditions for various gas mixtures. 


under wind-tunnel conditions are represented by a single 
curve since a fixed reservoir temperature (about 40° C ~ 100° 
F) has been assumed. A single mixture of 25-percent 
CBrFj-with 75-percent monatomic gas is considered in figure 
4 (a). These proportions are about right for M co = 2.0. At 
M m = 2.5, y a is 1.44 and y, is 1.34. This variation does not 
seem excessive. At M m = 1.3, y is uniformly less than 1.4, 
varying between 1.37 and 1.34. In view of previous com- 
ments about the apparent insensitivity of transonic flows 
to much larger variations in y, as well as the complete in- 
sensitivity of subsonic flows to variations in y, it would 
appear that a single mixture could be used satisfactorily 
from low subsonic speeds to at least = 2.5. 

The situation is more complicated within the Mach num- 
ber range between 2.5 and 3.5, as figure 4 (b) illustrates for 
the case of a single mixture of 40-percent CF 4 with 60-percent 
monatomic gas. It is to be noted that a variation of y is not 
important at conditions close to free-stream conditions. For 
example, y does not appear within the body of linearized 
subsonic or supersonic theory, but first appears in second- 
order terms. Hence, Busemann’s second-order theory for 
two-dimensional flow has been used to estimate the range 


within which y must be maintained in order to introduce less 
than 1 -percent error in A p/q. The boundaries of this range, 
indicated by dotted lines in figure 4 (b), shrink together as 
M„-» 1 and as M m - : > <» . For the moderate Mach numbers 
under consideration, it is seen that only the region of abscissa 
between about 0.2 and 1.0 is important when assessing the 
variations in specific-heat ratio. Consequently, it is de- 
duced that the gas mixture of figure 4 (b) approximates 
flight conditions at M m =3.5 about as well as does air in a 
wind tunnel. Such is not the case, though, at M a — 2.5 
(top portion of fig. 4 (b)), and it is not known whether 
moderate y variations near the stagnation region of the 
magnitude shown are important. 

If moderate variations in y are important, then the flow in 
an air wind tunnel will also show important differences from 
the flow in flight through air at Mach numbers of about 4 or 
greater. Throughout this report the operating temperature 
of a wind tunnel is assumed to be, as is current practice, 
either room temperature or that temperature which is 
sufficient to avoid condensation. Figure 5 compares the 
values of y in an air wind tunnel (long dashed curves) with 
the corresponding values of y in flight (cross-sectioned 
region). The observed y differences between air in flight 
and in wind tunnels are considerably greater than any of 
figure 4 (b) between flight and gas-mixture wind tunnels. It 

Gas-mixture wind tunnel, 0.4-CF4 + 0.6-Ar 

Gas-mixture wind tunnel, 0.2- CF 4 + 0.2-Ar + 0.6-Air 

Air wind tunnel 




o 


o 

cc 



Figure 5. — Variation of specific-heat ratio between free-stream and 
stagnation conditions at hypersonic speeds. 
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is also evident from figure 5 that the temperature dependence 
of 7 for polyatomic gases may actually represent a real virtue, 
particularly in hypersonic research. For polyatomic gases 
under wind-tunnel conditions, as well as for air under flight 
conditions, these caloric imperfections arise from a common 
physical phenomenon, namely, temperature dependent 
energy of vibration between atoms within a molecule. Under 
flight conditions at Mach numbers of about 6 or greater, 
additional energy can be transferred through dissociation, 
which is responsible for values of 7 less than 1.28 in flight. 
Clearly, the ratio of specific heats is one dimensionless pa- 
rameter pertinent to hypersonic flight through air — which is 
not simulated by air at wind-tunnel temperatures — but 
which is simulated by certain gas mixtures at wind-tunnel 
temperatures. 

MACH NUMBER LIMIT FOR CONDENSATION 

In order to achieve hypersonic Mach numbers in a wind 
tunnel without supersaturating the test medium, it is neces- 
sary, of course, to heat the medium. The approximate | 
reservoir temperatures (estimated to nearest 25° K) cor- j 
responding to the particular gas mixtures considered in 
figure 4, are as follows: 


Gas 

Reservoir temperature °K to 
prevent saturation; p t = 100 psia 


A/ro=5 

A'f co = 6 

Afro — 7 


350 

450 

550 

0.4 — CFt+O.G — Ar 

600 

675 

0.2— CF,+0.2 — Ar+0.6 air 

500 

625 

750 


As might be expected, the gas mixtures require higher operat- ; 
ing temperatures than does air. 

If is important to note that the trend of the results tabu- 
lated above can be reversed by not requiring that the 7 vari- 
ation of tire gas mixtures in a wind tunnel duplicate that of 
flight through air. Certain pure polyatomic gases require 
less heating to achieve hypersonic flow without saturation 
than does air, as the following table illustrates: 


Gas 

T„, °C 

Too 

7 t 

Reservoir tem- 
perature °K to 
prevent satura- 
tion at p, = 100 
psia; A/m = 7 

Air 

-180 

1.40 

1.38 

550 

CF< 

-128 

1.33 

1. 12 

450 

SFa 

— C>4 

1.27 

1.07 

450 

CF 3 CF 3 

-78 

1. 19 

1.07 

375 


For the three gases, air, CF 4 , and CF 3 CF 3 , the above trend is 
opposite to what would be expected from consideration 
solely of the respective boiling points (—180° C for air, 
— 128° C for CF 4 , and — 78° C for CF 3 CF 3 ) because of the 
dominating effect of the reduced values of 7. Reducing 7 
will reduce the rate at which temperature varies with pres- 
sure in an expansion process, since T~p (T_1)/T . This results 
in producing lower pressures, and hence higher Mach num- 
bers, before the temperature is reduced to a point where 
saturation occurs. From a physical viewpoint, this 
phenomenon can be traced to the conversion of internal 


vibrational energy of a molecule into directed kinetic energy 
during expansion of a polyatomic gas. To achieve a given 
ratio of directed to random energy (given Mach number), 
the polyatomic gas with large vibrational energy does not 
have to be expanded to as low a temperature as does a gas 
without vibrational energy. 

deviations from thermal equilibrium 

(HEAT CAPACITY LAG) 

When a polyatomic molecule passes through a region of 
rapid change in temperature, the energy stored in vibration 
between atoms within the molecules does not always adjust 
to its environment with sufficient rapidity to maintain 
thermodynamic equilibrium. This leads to a time lag in 
the heat capacity determined by the relaxation time r. 
Measurements show that r increases if either pressure or 
temperature decreases. An excessively long relaxation time 
could significantly affect a high-speed flow, especially in a 
wind tunnel where the pressures and temperatures encoun- 
tered are low. 

Detailed computations are presented in Appendix B which 
provide an estimate of relaxation-time effects for various 
gases. These estimates are based on a relation between 
relaxation time and bulk viscosity. The derivation of this 
relation, and a discussion of its range of validity are presented 
in Appendix C. Here, only end results are discussed. The 
relaxation time (r a ) at atmosphere conditions is known for 
some of the fluorochemical gases, and generally is the order 
of 10 -7 second. For others, r a is not known from direct 
measurements, but can be estimated from related measure- 
ments by using certain empirical relationships between 
relaxation time and molecular structure. Thus, t„ for 
CBrF 3 is estimated to be in the range 1CT 8 to 10~ 7 second, 
and that for CF 3 CF 3 to be in the range 10~ 9 to 10“ 8 second. 

Relaxation time is not the only important quantity which 
must be considered. The basic parameters which measure 
effects of heat-capacity lag on a gas flow involve the product 
tC u where (\ is the heat capacity of internal vibrations 
within the molecules. From the estimates in Appendix B, 
heat-capacity lag appears of more importance in altering 
boundary-layer flow than in altering pressure distribution, 
and is of greatest importance at low Reynolds numbers. For 
air under wind-tunnel conditions, r is so long that the small 
internal vibrational energy effectively is “frozen,” thereby 
rendering relaxation effects negligible even at low Reynolds 
numbers. For most of the gas mixtures under wind-tunnel 
conditions, the estimated relaxation effects arc negligible at 
Reynolds numbers of the order of 10 5 or greater. They can 
be significant, however, at low Reynolds numbers (the order 
of 10 4 or less). It should be remembered that for air under 
flight conditions, C t is not negligible as it is in a wind tunnel. 

’ Thus, relaxation effects can be significant in flight, although 
the estimates of Appendix B suggest this may be the case 
i only at high altitudes (low Reynolds numbers) and at high 
flight speeds. In broad terms, heat-capacity lag is esti- 
mated to be significant at Reynolds numbers below about 
10 4 for moderate supersonic speeds, and below about J O 6 for 
hypersonic speeds. Although relaxation-time effects in 
flight would not be simulated in an air wind tunnel, it is 
possible to simulate them approximately in a gas-mixture 
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wind tunnel by employing two polyatomic gases in the 
mixture — one having rCi greater than for flight through air, 
and the other having tC { less than for flight through air 
This possibility appears to be of interest in low-density, high 
Mach number research, and provides another example of an 
effect pertinent to hypersonic flight through air which cannot 
be studied in an air wind tunnel but can be studied in a gas- 
mixture wind tunnel. 

It should be mentioned that the idea of employing more 
than one polyatomic gas in a mixture might be utilized profit- 
ably in other ways. For example, various pairs of some poly- 
atomic gases when mixed in special proportions form an 
azeotropic mixture; that is, a mixture for which the boiling 
temperature is not intermediate between the boiling tem- 
peratures of the two constituent gases, as ordinarily is the 
case, but is lower than the boiling temperature of either con- 
stituent. A mixture of 32-percent C 3 F S with 68-percent 
CH 2 F 2 , for example, boils at —58° C, whereas the respective 
individual gases boil at —38° C and —52° C (ref. 63). Thus, 
although neither C 3 F 8 nor CH 2 F 2 could individually be used 
at M«=2.5 under normal wind-tunnel operating conditions 
without danger of condensation (and, hence, have not been 
included in table IV (b) ) , their azeotropic mixture would be 
usable under such conditions. Other gas pairs known to 
form azeotropic mixtures are CHF 3 with C 2 F 6 , CFoCL with 
CF 2 CFCF 3 , CHC1F 2 with CF 2 CFCF 3 , and CF S C1, with 
CH 3 CHF 2 (see refs. 48 and 63). A related point worth 
mentioning is that certain polyatomic gases which would not 
be usable singly because of flammability characteristics may 
be usable when mixed with other gases. An extreme example 
illustrating this point is the obviously flammable butane 
(C 4 H 10 ) which, when mixed with CF_C1 2 in portions up to 
about 30 percent, is no longer flammable (ref. 64). 

VISCOSITY VARIATION WITH TEMPERATURE AND PRANDTL NUMBER FOR 
VARIOUS GASES 

Thus far, mention has been made of two effects associated 
with hypersonic flight through air which are not simulated 
in an air wind tunnel operating at conventional temperatures 
but which can be simulated in a gas mixture wind tunnel. 
These effects pertain to specific-heat variation with tempera- 
ture and to relaxation-time phenomena (or bulk-viscosity 
phenomena). To this can be added a third effect; namely, 
that pertaining to variation of viscosity with temperature. 
It is to be remembered, though, that not all of the three 
effects are significant simultaneously in a given research 
problem; likewise, not all can be simulated simultaneously 
by a single gas mixture. One gas which — in a wind tunnel — - 
approximates the temperature-viscosity relationship of 
flight conditions is pure helium, as illustrated in figure 6. 
A different gas which approximates the corresponding rela- 
tionship of air in a wind tunnel is pure argon, as also illus- 
trated in figure 6. The various curves representing wind- 
tunnel conditions in this figure correspond to reservoir tem- 
peratures which will avoid saturation at M= 5. Helium and 
argon have the same value of y, and very nearly the same 
value of Pr, they differ in aerodynamic behavior only in 
then 1 viscosity-temperature relationship. Separate measure- 
ments with these two gases, therefore, would single out any 
effect of viscosity-temperature relationship on a given 



M m —5. 

phenomenon. This possibility would appear to be useful in 
boundary-layer research. 

When gases are mixed in arbitrary proportions, the Prandtl 
number can vary markedly. A rather extreme example of 
this is afforded by He and SF 6 mixtures. At 0° C the 
Prandtl number of He is the same as that of SF 6 , namely, 
0.69. Mixture of these two gases, however, can result in 
values of Pr as low as 0.21, according to the calculated curve 
shown in figure 7. Similarly, Pr for argon is 0.67, but mix- 
tures with helium can result in values as low as 0.41. These 
mixtures provide values of Pr less than that of air (0.7 at 
room temperature), but other mixtures can provide values 
greater than air. For example, mixtures of CC1 2 F 2 and He 
(results not shown) cover the Pr range between 0.82 and 0.22. 
As a test of the calculation method, which indicates such large 
variations in Pr, some experiments are available for the case 
of He-A mixtures (ref. 65 for viscosity, ref. 66 for conductiv- 
ity). Calculation and measurement agree well for this case 
as is evident from figure 7. The particular equations em- 
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ployed in the calculations 5 are the rigorous mixture equa- 
tions of kinetic theory developed by Hirschfelder, et al. (ref. 
67). It is noted that the large reductions in Pr can be attrib- 
uted primarily to the mixing of light with heavy molecules 
(m 2 /mi = 37 for SF 6 -He, mojirix — 10 for A-He). If gases of 
more comparable molecular weight were mixed, say A and 
CF. t (m a /m 1 =2.2), the corresponding reduction in Pr would 
not be large. Large variatious in Pr afforded by special 
mixtures of gases would appear to be a significant tool for 
conducting supersonic heat-transfer research, since the 
Prandtl number is an important parameter in such phe- 
nomena. 

CONCLUDING REMARKS 

By employing in place of air a test medium consisting of a 
heavy monatomic gas mixed with a heavy polyatomic gas in 
proportions that yield the desired value of 1.4 for the ratio 
of specific heats, it appears possible in conducting aero- 
dynamic research to realize certain significant advantages. 
The existence of some of these advantages lias been demon- 
strated by previous investigations concerned with pure 
heavy, polyatomic gases, especially the Freons. Because 
pure polyatomic gases have a specific-heat ratio much 
lower than that of air, they cannot be employed to yield 
aerodynamic data directly applicable to air for flows involv- 
ing extensive supersonic regions. Aerodynamic data ob- 
tained with the gas mixtures considered herein would 
directly apply to air. 

Perhaps the most significant result of this study is the 
observation that a gas mixture can be concocted which 
behaves — under low-temperature wind-tunnel conditions — 
dynamically similar in several respects to air under high- 
temperature flight conditions, whereas air under low-tem- 
perature wind-tunnel conditions behaves dissimilarly. This 
situation arises because the flow of two gases can be made 
dynamically similar on a macroscopic scale, even though 
they differ in microscopic structure, provided all the per- 
tinent, dimensionless, macroscopic parameters (such as 
involve, for example, relaxation-time phenomena, tempera- 
ture variation of specific heat, and temperature variation of 
viscosity), are duplicated between the two gases. In 
achieving dynamic similarity, the dimensional quantities, 
such as temperature, are not duplicated. 

Many of the fluorochemical gases studied in the present 
research are more inert than Freon 12. This extreme 
inertness is highly desirable, especially since some of these 
gases when mixed with the necessary small amount of 
oxygen can be breathed even in large quantities without 
observable effect by animals, and presumably also by man. 

5 These calculations were made by Mr. Donald Doran. Intermoleeular forces were rep- 
resented, in the nomenclature of reference. 67, by a Lennard-Jones potential for He-SFe and a 
modified Buckingham potential (with a = 12.4) for Hc-A. Force constants e/lv were 9.16“ K 
for He, 239° lv for SFa, and 123.2° K for A. Molecular dimensions were selected to yield the 
proper values of each transport property of the pure gas component. 



Figure 7. — l’randtl number variation for two gas mixtures containing 
helium; T= 0° C. 

Likewise, because of their inertness they could be used 
without undesirable effect on machinery. 

Since the magnitude of the advantages that can be 
realized through use of gas mixtures has been shown to be 
large, they may outweigh the evident operational dis- 
advantages which result from use of any gas other than ah’. 
Whether or not the advantages surmount the disadvantages 
can only be ascertained by detailed study of individual 
cases. Clearly, though, the possible use of gas mixtures 
warrants consideration in the design of new research facilities. 
In view of the magnitude of the advantages, it appears that 
some experimental work with gas mixtures is in order. 


Ames Aeronautical Laboratory 
National Advisory Committee for Aeronautics 
Moffett Field, Calif., Mar. 12, 1954 



APPENDIX A 


ESTIMATE OF DEVIATIONS FROM THERMAL PERFECTION 

By definition, a thermally perfect gas obeys the equation 
of state pm—pRT given previously as equation (2). Real 
gases closely follow this equation if the density is low, or 
if the temperature is near the Bojde temperature (about 
2.6 T e ), but deviate from it under other conditions. A 
more exact state equation is 

5 l f =l+ RT B{T) (A1) 

where B(T) is the second virial coefficient having dimensions 
of volume. According to the law of corresponding states, 
B(T) is proportional to the critical volume V c . 

From the results of reference 67, the second virial coefficient 
for a pure gas is approximately 0.75 V C B*(T*), where T* = 
\.3TIT C and where B*(T*) is a fixed dimensionless function 
tabulated in reference 67. For a mixture of two gases the 
appropriate relations are : 

5(7 7 )=x 1 2 5 I (T 1 )+2r 1 a: 2 fi 12 (2 7 )+r 2 2 5 2 (r) (A2) 

where 

B-0.75V cl B*(1.3T/T cl ) 
B 2 =0.75V c2 B*(l.3T/T c2 ) I (A3) 

/T7 1/3 I T/ l/3\ 3 

B n =0.75( J cl y c2 ) B*(] .3T/-\ T C] T c2 ) 

By use of these equations, the term pB(T)/RT of equation 
(Al)— which represents the fractional deviation from 
thermal perfection — has been calculated for the gas mixtures 
listed in table IV. Computations were made both for 
reservoir conditions ( p t ,T ,) and test-section conditions 
( p m ,T m ) at Mach numbers of 1.3, 2.5, and 3.5. For M„ = 
1.3 the values of p t B{T)IRT t and p a B(T a )jRT a were about 
the same. The average value of \pB(T)/RT\ per atmosphere 


of reservoir pressure varied from about 0.0014 (C 4 F 8 , C 3 F 8 , 
CF 3 CF 3 , CF 4 , and NF 3 mixtures) to about 0.004 (CH 2 F 2 mix- 
ture). Inasmuch as transonic research apparatus commonly 
are limited to reservoir pressures of several atmospheres (often 
because of critical model loads), the resulting deviations 
from thermal perfection would be within about 1 percent 
and can be neglected. For M m = 2.5 and M„ = 3.5, the 
values of | />„,/?( 7^ )/RTJ are substantially lower than 
\p,B(T,)/RT[\. Both are tabulated as follows: 


1 

' .Vf„ 

=2.5 

i 

i JVf CO 

=3.5 

Gas 

l p oo b ( r co ) I 

\P.B(T,)\ 

|pa>B( Tco)| ! 

\P,B(T,) | 


RT m 

1 

RTi 

Too 

RTi 1 


CBrFs... 0.0017 0.0033 

SFr .0013 . 0027 

CF.CF. . .. .0010 . 0018 

CCIFi... 0010 . 0030 

CF, 0012 . 0019 0.0007 0.C023 

N’Fj 0012 . 0019 . 0007 . 0022 

CHFj .0017 . 0034 

Air ... .0005 . 0002 . 0002 . 0002 


It may be deduced that gas mixtures at test-section condi- 
tions (p„,T„) would deviate 1 percent from thermal 
perfection at reservoir pressures between about 6 and 10 
atmospheres for M a = 2.5, and at about 14 atmospheres for 
M m = 3.5. Since practical operation pressures of wind 
tunnels are within these limits, the flow of these gas mixtures 
over a typical model can be regarded as closely approximat- 
ing that of a thermally perfect gas. Deviations from thermal 
perfection at reservoir conditions ( p t ,T ,) are two to three 
times as large as at test-section conditions. Hence, when 
computing test-section Mach number or dynamic pressure 
from measured values of p t and T,, the deviations near 
reservoir conditions may have to be considered for reservoir 
pressures of the order of 5 to 10 atmospheres or greater. 
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APPENDIX B 


RELAXATION-TIME DATA AND ESTIMATE OF EFFECT FOR 
VARIOUS GAS MIXTURES 

Results of measurements of relaxation time at atmospheric 
pressure (t„) for some of the polyatomic gases considered in 
this report are as follows: 


Oils 

Temper- 

ature, 

°K 

To 

micro- 

seconds 

Refer- 

ence 

CCI 2 F 2 

293 

0. 09 

68,4 

! CF, 

373 

.00 

69 

SFs 

293 

.59 

70 

CHFj. 

373 

.42 

09 

CHCIF 2 .. 

293 

.10 

68 

CH 2 F 2 

373 

.055 

69 

CO" 

293 

7.0 

08 

N 2 0 

293 

1.1 

68 


Although measurements of r a could be found only for those 
gases listed, certain known characteristics of relaxation time 
enable a rough estimate to be made for some of the remaining 
gases; for example, molecules with freedom of internal 
rotation invariably have very short relaxation times, gen- 
erally less than 10~ 8 second (see ref. 71, for example). Also 
the smaller the lowest fundamental vibration frequency, the 
shorter the relaxation time (see refs. 68 and 69). These 
general characteristics indicate that fluorocarbon molecules 
with internal rotation (e. g., CF 3 CBrF 2 , CF 3 CC1F 2 , CF 3 CF 3 , 
CH 3 CC1F 2 CH 3 CHF 2 , C 3 F s , and C 4 F 10 , but not C 4 F 8 ) would 
have relaxation times (r a ) in the range 10 -8 to 10 -9 second, 
and that the methane-like gases CBrF 3 and CBrClF 2 
would have a relaxation time in the range of 10~ 7 to 10~ 8 
second. It is to be noted that relaxation times listed are 
for pure gases, and that the presence of certain impurities 
can reduce r a for those gases having relatively long relaxation 
times. 

Inasmuch as the relaxation time depends on temperature 
and pressure, some method of estimating r for the par- 
ticular conditions encountered in wind tunnels must be 
employed. The analysis of Bethc and Teller (ref. 72) yields 
for the relaxation time r„ of the lowest frequency (/) mode 


Z 10 

T ' Z(l~e~ h ' /kT ) 


(Bl) 


where Z is the total number of collisions one molecule 
experiences per second, and Z l0 is the average number of 
collisions required to deactivate the lowest mode from the 
first quantum state to the zero state. Subsequent calcula- 
tions are based on the additional equations 


Z=1.3 - 


(B2) 


which is a result from the rigorous kinetic theory of gases 
(ref. 73), 


r 


CjT , 

c. 


(B3) 


which is a result of the assumption (ref. 74) that once the 
lowest frequency mode (specific heat C r ) is excited, energy 
is then rapidly transferred internally to all other modes 
(total internal specific heat C\), and 

Z 10 ~ 25(~)-' ! e b/TM (B4) 

which is an approximate result of Bethe and Teller (ref. 72). 
In this last equation, 6 is a constant to be evaluated for 
each gas from measurements of r a at the temperatures listed 
in the above table. These equations were developed for 
pure polyatomic gases. They 'will be used also for gas mix- 
tures since the basic theory (especially eq. (B4)) is highly 
approximate, and since the effectiveness of an argon mole- 
cule in exciting vibration generally is roughly the same as 
that of an inert polyatomic molecule. (See the summary 
table in ref. 75 and the specific data for CHC1F 2 -Ar mix- 
tures in ref. 76.) 

Relaxation times r„ corresponding to free-stream condi- 
tions (TnjPn) in a wind tunnel have been calculated for 
assumed reservoir conditions of one atmosphere total pres- 
sure and 40° C total temperature. The results, together 
witli the “relaxation distance” for the gas mixture 

were computed to two significant figures and are tabulated 
to one significant figure as follows: 6 


Polyatomic 

component 

M ® =1.3, 
Iv., pee = 

Too =230° 
0.30 atm 

Moo =2.5, Tco =140° 
K, p 00 =0.059 atm 

Mas =3.5, Tco =105° 
K, poo =0.013 atm 

r co , sec 

d, in. 

. 

To,, sec 

d , in. 

Tea, sec 

d t in. 

CCI 2 F 0 

2X10-7 

0. 002 





CF, 

2X!0-« 

.02 

2X10-’ 

0.3 

X 

0 

3 

SFa 

1X10-° 

.01 

9X10 -6 

.15 



1 CHFj 

2X10-° 

.02 

1X10- 3 

.2 



CHCIF 2 

2X10-7 

.002 





CH 2 F 2 

9X10- S 

.001 

T 

0 

X 

*0 

.008 



1 C0 2 

3X10-5 

. 4 




, N 2 O 

4X10-0 

.05 


1 


i 




l 



It is to be remembered that under wind-tunnel conditions, 
air has a much longer relaxation time than any of the gases 
tabulated above. At room temperature, for example, r for 
air is the order of 10~ 3 second (ref. 77), which is 10 2 to 10 4 
times as long as for the polyatomic gases. This does not 
mean, however, that relaxation effects are important for air 
in a wind tunnel because the amount of internal vibration 
heat capacity ((?<) that lags is very small. The effect of 
heat-capacity lag on flow conditions is influenced by both 
Cj and r. 

To determine the dimensionless parameters which form a 
pertinent measure of relaxation effects, use is made of the 
theoretical equivalence between relaxation phenomena and 

0 Some of these values differ from corresponding values tabulated in the original publication 
of this research (TN 3226) . Such differences arc not large and are due to the use of a different 
definition of relaxation time (t) in equation (Bl). The relaxation time originally used is 
not consistent with the definition of Bethe and Teller, whereas that used above is (see 
Appendix C). 
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volume-viscosity (k) phenomena. This equivalence is do- 
rived and discussed in Appendix C where the fundamental 
relation K=})T('Y—l)C i /C v is shown to be valid under condi- 
tions where the relaxation time is small compared to the 
characteristic time (l/u^ for aerodynamic problems). For 
purposes of estimation, then, the complete, first-order stress 
tensor, including the volume viscosity, is used as the basis 
of computation. 

PiMVnh„, Kr +^ a &u (B5) 

•Stoke* ** « 

In this equation k is the coefficient, of volume viscosity, 
du a /dx a is (using summation convention) the divergence of 
the velocity vector, and is unity if i is equal to j, but is 
zero otherwise. The above equation, together with the 
equation relating k to r, will be employed to estimate effects 
of heat-capacity lag both for inviseid flow and viscous flow. 

For inviseid flow the Xavier-Stokes portion of the stress 
tensor is simply the static pressure p. The pertinent simi- 
larity parameter involving relaxation time is nu x /(pl) — 
Since the equation of steady motion 
involves pUjdu,/dxj and bp,j/bxj, the ratio of the relaxation 
term to the other stress term is, in order of magnitude, 



which is seen to be the same as the similarity parameter . 7 
Hence, nu^/pl is interpreted as a measure of the fractional 
effect of relaxation phenomena on pressure distribution. 
In a gas mixture, Computations of the resulting 

parameter Kiu/pJ=y(y—l)(T a> u <c /l)C,/(C l> ) for wind- 
tunnel flow conditions over a 1 -inch model yield the following 
results: s 


r>C x 2 T<n Hen _K» co for /= 1 ill. Sllld 


Polyatomic ^ $ yj/- 1 atm 

component — 

i 

A/co = 1.3 A/ oo = 2.5 1 A/co =3.5 


rcizixj . o. ono2 ( 

ri'i .002 ' 0 01 0.04 

SIV .002 .013 

nil', .001 . oof! 

, CHC1F_> .002 

' .00005 .00008 

i oo 2 .02 

XiO .003 


For transonic wind tunnels it is seen that only C0 2 mixtures 
would exhibit appreciable relaxation effects on pressure 
distribution. On the other hand, for .\/„ = 3.5, mixtures 
involving CF., would he affected roughly by 4 percent for 
the assumed conditions of p, = 1 atmosphere and /=1 inch. 
To maintain these effects below about 1 percent would 
require the product p,l to be greater than 4 atmosphere- 
inches (e. g., a 1 -inch model in a 4-atmosphere tunnel). 

1 A more precise calculation of this ratio of stress linns can be made using tile continuity 
(■(ination to eliminate oiiVOr*. This pioccdurc ultimately yields a value 0cUa,/p/)|.\/3o(«/ 
" co)M r ll)\, from which it can he deduced that, the simple parameter («.«<„//)/) will tmder- 
c-dhnnlc relaxation elTccts for blulT obstacles like tire sphere, but will overestimate them for 
slender obstacles like an airfoil. 

s The contribution of rotational energy to * is disiegardcd since this is small compared to 
the contribution of vibrational energy. 


Such restrictions arc believed not to represent significant 
limitations on the usefulness of gas mixtures, especially 
since the favorable effect of certain impurities in reducing - 
has not been included in the analysis. It is concluded, then . 
that insofar as pressure distribution in inviseid flow is 
concerned, relaxation effects need not be considered for the 
usual wind-tunnel operation conditions at Mach numbers 
near 3.5 or below. 

For viscous How within a laminar boundary layer, the 
Xavier-Stokes portion of the stress tensor customarily is 

^ XL 

approximated by p The ratio of the relaxation term to 
the ordinary term in the equation of motion is 

b / bu k \ 1 / bu a \ bu a 

dr \ dr,. / / \ bx a ) / 5 \ k d ,r a 

d / dw.\ 1 / d u\ \J.) p bn 

by\b y) i y d y) d ij 

Examination of the solution to boundary-layer equations for 
flow over a flat plate (e. g., ref. 78) indicates the average value 
of (bu, : /bx k )/(biijby) to be of the order of (y — 1 ) -l /,* 2 1 0 A lit . 
and 6/7 to he of the order of [G-f (7— 1 ).)/„■’] ylh. Hence. 


d / 

{ . du ") 


d.r ' 

\ b.r a ) k 

-(7-1).1/ x H 

5 I 

( bxi\ p 

10 /A J 

by 

K^by) 



■J 10+ Cr- 


1 )-!/,. "1 


(B7) 


This ratio, which is interpreted as a measure of the fractional 
effect of relaxation phenomena on boundary-layer flow, is 
important primarily at low Reynolds numbers and high Mach 
numbers, just as is the effect of relaxation phenomena on 
inviseid flow. 

From the equation KU^pl — Up) (yM^-Re) and from the 
preceding table, the ration k'p corresponding to froe-stream 
conditions can he computed for the various gas mixtures. 
The results are as follows: 


Calculated values of (a » 

Polyatomic com j >oncnt 

A/ x = 1.3 A/ * =2.3 A/co = 3.5 


CC1>F_. 37 1 

C F i _ _ 350 410 370 

$ F* 300 3S0 

CI1 F-> 230 170 

CIIC1K* 37 ... 

Cn>F*. - S 2 

CO» 3700 

N >0 530 


At stagnation conditions njp would be the same order as that 
tabulated for d/ a> =1.3; hence, broadly speaking, the various 
fluoroehemieal gas mixtures cover a range in values of k p 
between about 2 and about 400. For the gases of most 
interest, namely, CF 4 and SF 0 , the corresponding value of 
k ! p is seen to be about 400. Consequently, the value n'p— 
400 can be substituted into equation (B7) to obtain an 
estimate of relaxation-time effects on boundary-layer flow. 
The result is that the beat-capacity lag is found not to be 
significant except at low Reynolds numbers. In particular, 
at a Mach number of 3.5, it is indicated that 011 I 3 - below or 
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near Reynolds numbers of about 2,000 would the effects of 
heat-capacity lag be significant. Even at a Mach number of 
7 they would be significant only near or below Reynolds num- 
bers of about 20,000. 

Turning now to the case of air, a distinction must be made 
between wind-tunnel conditions and the conditions of 
flight.. In a wind tunnel, the relaxation time for vibration 
is so long that the essential assumption inherent in the 
k—t relationship is not a valid assumption, since the vibra- 
tion energy can lag far behind the translational energy. 
The energy in vibration for air at wind-tunnel temperatures, 
however, is negligible so that there are no appreciable 
effects of lieat-capacitv lag under such conditions. In 
flight, the situation is different, since the higher tempera- 
tures encountered reduce the relaxation time to values 
which are low enough to enable the vibrational energy 
almost to follow the translational energy. At a tempera- 
ture of 1500° Iv, for example, the experiments of Blackman 


(ref. 77) show the relaxation time for air to be the order of 
10 -5 second. At 3000° K it is the order of 10 -0 second. 
These times arc short enough to apply equation (C9) at the 
temperatures of flight. The result is that the ratio x/ti is 
calculated to be the order 1,000 at 1500° K and the order 
of 50 at 3000° Iv. Even if the higher value (k/h = 1000) is 
used and substituted into equation (B7), the effects of heat- 
capacity lag in flight are found to be small, except under the 
combined conditions of a high Mach number and simultane- 
ously a low Reynolds number. Under such conditions, the 
heat-capacity lag could be significant. In order to simulate 
in a low-density gas-mixture wind tunnel the licat-capacity 
lag that would exist in the boundary layer during flight, 
the ratio kIii would have to be simulated. It is possible to 
simulate approximately this ratio in a gas-mixture wind 
tunnel inasmuch as the range of values of k/h covered by 
the various polyatomic gases bracket the corresponding 
values for flight conditions in air. 



APPENDIX C 


RELATION BETWEEN BULK VISCOSITY AND RELAXATION 

TIME 


Tito calculations of Appendix B arc based on a relation 
between relaxation time r and bulk viscosity coefficient k. 
This relation can be derived from the general equations of 
motion for a viscous compressible gas by considering the defi- 
nition of relaxation time and one of the fundamental prop- 
erties of entropy. The method of derivation presented here 
is rather general and leads to a simple result which is com- 
pared later with results of other analyses of the problem. 

The general equations of motion involve, the viscous stress 
tensor 


/but buj 2 bu„\ bu a 

ysj; ,+& 7+3 5 « 57. ; +t 5 7 s <- 


~(pu\. ■ 

\ / Nrtvicr 

Stokes 



(Cl) 


and the energy dissipation function 


_1hj /d«j dttA 

2Aa5 ; +57,) 


2 \dxj ' ds t 

a >««y 


= *P X :i v ier -4- K ( 

Stokes 




= ^X.’ivier-J - 
Stokes’ 


\dx 


(C2) 


where the operator T) represents the substantial differential, 
and the Xavier-Stokes portion of ^involves only the ordinary 
coefficient of viscosity. The bulk viscosity k — just like the 
ordinary viscosity p — is regarded as a state property and, 
therefore, can be evaluated as a function of temperature and 
pressure from any pertinent experiment. For example, 
suppose an experiment with Court te-type flow were selected 
wherein the effects of bulk viscosity are negligible and the 
effects of ordinary viscosity are large (as would be the case 
at low speed under which condition dv a fd.r a = — Dp/pDt is 
negligible). The viscous stress tensor would reduce to 
p rj — pdu/dy, so (hat measurement of the shear and the 
velocity gradient would enable the state property p to be 
determined. On the other hand, suppose an experiment 
were selected wherein (he effects of ordinary viscosity are 
negligible and the effects of bulk viscosity are large (as 
would be the case for the rapid compression in front of a 
small pitot tube or for the velocity dispersion of ultrasonic 
waves, under which conditions the time rate of change of 
density is the dominant term). In this case, the dissipation 
function (eq. (C2)) simplifies to the single term n(Dp/pDt ) 2 , 
so that measurements of the energy dissipation (total pres- 
sure loss) and the rate of density variation would enable the 
stale properly k to be determined. 

By analysis of a flow wherein DpjDl is unusually large and, 
hence, wherein the bulk viscosity term dominates all other 


terms in the stress tensor, a relationship can be established 
between k and r. For (his type of flow, the energy equation 
can be written 


P T 


DS k ( D P \ 
Dt ~ <f —p' 1 \Dt ) 


(C3) 


where S is the irrerciviblc entropy increase. The definition 
of relaxation time r is 


DE t _E,(T)-E t 

Dt T 


(C4) 


where T is the temperature corresponding to the transla- 
tional degrees of freedom, E^T) is the equilibrium value of 
internal energy corresponding to the temperature T, and E t 
is the actual internal energy (corresponding to equilibrium 
at a different temperature 77) . From thermodynamic theory 
a third equation is obtained which completes the building 
blocks of the analysis; it provides an equation for the change 
of entropy when an amount of energy E t flows from one 
source (translational energy at temperature T) to another 
source (vibration energy at temperature T t ). 


DS= 



JDEj 

T 


The above three basic equations can be combined to yield 
the desired relation between bulk viscosity and relaxation 
time. 

It is emphasized that no assumption is made about the 
variation of specific heat with temperature, nor is any as- 
sumption made either about the magnitude of the internal 
energy or about the form which that internal energy takes. 
For example, the internal energy could be in rotation of 
the molecules, or in vibration between the atoms of molecules, 
or in electronic excitation. The essential assumption made 
is that the internal energy E, does not lag far behind the 
equilibrium internal energy corresponding to the tempera- 
ture T. In other words, it is assumed that the relaxation- 
lime' effects arc small, so that a first-order theory is appli- 
cable. Tlie various quantities are illustrated in figure S. 
Since a first-order theory is being considered, the internal 
specific heat C t can be regarded as constant over the small 
temperature range 7 1 — 7’,, and the following approximation 
can be made: 


DS—DEt {T Ei( ^ Ti E< 


(C6) 


By combining this entropy equation with the equation de- 
fining relaxation time, there results 

DS_mj)-E t Y 
Dt tT-C\ ~ 
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Because the special experiment selected for consideration is 
ail isentropic one except for first-order irreversible processes 
associated with k, the following; additional equations result 
from equation (C:i ) : 

pT DS_ 1 / Dj> y 1 (I)'J'y 

T DC' P - \lJt ) ~(7-l fP V Di I 


— __ - c -T- 
(7-11* ' 


DI: 
\ Dt 


*(: 


?.y_r 

t ) L(7-l )(\Tt_ 


(C'Sj 


The final relation between bulk viscosity and relaxation time 
is obtained by combining equations (CS) and (C7) (and 
noting that p = pli T—p(y— 1 )C V T) 


K =p T (t 1) ~v~ (C9) 

* V 

Inasmuch as C\ can represent either rotational energy or 
vibrational energy, several comparisons can be made with 
other analyses. In reference 79, Tisza compared the ap- 
proximate equations for the dispersion of ultrasonic waves 
with corresponding equations from the. Xavier-Stokes rela- 
tions and obtained exactly equation (C9). However, in 
comparing with the approximate equations from ultrasonic 
theory, Tisza assumed in elfect that the amount of energy 
in internal motion E t is small compared to the total energy, 
that is, he assumed that C\ is small compared to C r . The 
above analysis shows that l lie relation obtained by Tisza 
actually applies to the case where the internal energy is a 
large portion of the total energy. Tisza assumed, of course, 
that the departures from equilibrium are small enough to be 
represented by a first-order theory (this assumption is made 
in all of the analyses discussed). 

A different approach to the problem is taken by Kohler, 
reference 73, who employed the equations of kinetic theory 
involving the Maxwell-Bolt zmann equation. Kohler con- 
sidered two cases — first, that of rotational internal energy, 
and second, that of vibrational internal energy. The 
definition of relaxation time for rotation (r r ) employed by 


Kohler, however, is not the same as that employed above. 
The two are related by the equation 

C\ + Tr 

/F T 7 

v / 

where C t =3li/2 is the specific heat due to translational 
degrees of freedom. Letting / r be the number of rotational 
degrees of freedom, then C r — L/?/ 2, and Ii/(y — 1 ) = C r — 
(3 -f-/,) R/2, so that equation (C9) becomes 


which is identical to an equation developed by Kohler. 
When Kohler considered vibrational energy, ho used the 
same definition of relaxation time ns that used herein, but 
he assumed that the internal energy in vibration ( E rib ) is 
small even compared to the rotational internal energy. 
With this approximation Kohler could write C,+ C ri( ,« C r , 
and thus obtained the equation 


k = (7— 1 )]>t ?ryf>V 

r ^ rib 

Since (\—C C r within the approximation made by 
Kohler, his result for vibrational internal energy is consistent 
with the basic equation (C9) of this appendix. 

In reference 80 Meixner considers the problem from an 
entirely different viewpoint, namely, that of the theory of 
irreversible thermodynamics. Meixner assumes in effect, 
as did Kohler, that the internal energy is small compared to 
the translational energy. The equation Meixner obtained 
can be reduced to 

C t 

K~pr(y 1 ) p, q 

Just as in the case of Kohler’s analysis, this result is seen to 
be consistent with that of equation (C9), provided it is 
remembered that the assumption of small internal energy 
enables C t ~Ci to be approximated by C r within the frame- 
work of Meixner ’s analysis. 

In reference 81 Wang-Chang and Uhlenbeek considered 
the problem from the viewpoint of kinetic theory, only 
they did not assume as Kohler did that the. internal energy 
was small. They did assume, however, that the internal 
energy did not vary significantly with temperature; that is, 
they assumed that the internal specific heat, was constant, 
as th is simplified their equations. Their final result can be 
written in (he form 

K =p T ( 7 — 1 ) ~ 

v r 

which is identical to equation (C9) above. Thus, wc see 
that the end result of the analysis of Wang-Chang and 
Uhlenbeek, just as the analysis of Kohler for rotational 
internal energy, actually applies to more general conditions 
than they assumed. Consequently, the derivation given 
above is seen to be rather general and to be confirmed by 
four separate analyses employing methods from several dif- 
ferent fields of science. 
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TABLE I. — PROPERTIES OF HEAVY .MONATOMIC GASES AND AIR CONSTITUENTS 


[ 



Boiling 

Critical constants 

Viscosit v 

Symbol 

Oas 

m,- r , • 
mole 

temp., 

° C 

T„ » K 

Pc. atm. 

l cc/gin. 
mole 

at 0° C ti 0 . 
micro poise 

Ar 

argon 

•JO 

-ISO 

151 

4S 

75. 3 

210 

Kr 

krypton 

S4 

-153 

210 

54 

92. 1 

230 

Xe 

xenon 

131 

-107 

200 

oS. 2 

114 

210 

x> .. 

nitrogen 

2S 

— 190 

120 

33. 3 

90.0 

1GG 

O; 

oxygen 

32 

-183 

154 

49. 7 

74.4 

189 


air 

29 


132 

37.2 

S3 

172 
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TABLE II.— GENERAL CHARACTERISTICS OF POLYATOMIC GASES 


Formula 

7712, 

gm. 

mole 

Chemical name 

Common or 
trade name 

Gcueral 
class ifi- 
cation 
( 1 ) 

Kcm arks 

References 

Spectro- 

scopic 

data 

Toxicity 

data 

Physical data 


238 

n-nerfluorobutane _ _ 

s 



14 



221 



S? 



50 

14, 50,44, 46,23 






Chemistry resembles NF 3 . Stable to at 







least 400° C. 




C 4 Fs 

200 


Freon C318 

s 


45 


48,15,26 








CFjCBrFj 

199 

bromoethforane 


s? 





CF 3 SF 5 

196 

trifluoromethylsulfurpentafluoride . _ 


s? 


— 



SeFi 

193 

selenium hexafluoride. 


s? 


50 


14,25,42 









188 

perfluoropropane ... 


s 



14 


C 2 F 8 NF 

171 

perfluorodimethylaminc . 


s? 



44,47 ] 

I 






known. 



CBrCIFa 

105 

bromochlorodifluoromcthane 

Freon 12B1. ... 

q 

Toxicity probably excessive. More toxic than 

49 

21 

2G 






CHClFj. 



i 

CFjCCIFj 

155 

chloroethforane . _ 


s 


33 



CFjOCFj 

154 

perfluoroether... . . 


s? 



50 

14, 56, 44, 23 






Liquid density unknown. Stable to at least 







400° C. 



! 

CBrFj. 

149 

bromotrifluoromothano 


s 


35 

21,56 

14, 2G, 56 




Kulene 131. 




SF, 

146 

sulfurhexafluoride 

1 

S, I 

: Excellent high-voltage gaseous insulator. Ex- 
tremely inert chemically. Thermally stable 

29 

17 

14, 56, 25, 27, 
' 52, 53 





1 


[ 


to about 500° C. 




' OF 3 CF 3 

138 

ethforane 

i 

1 s 

j Thermally stable to about 700° C_. 

i 34 


14, 5G, 19, 39, 22 
14, 26 

CCIF 2 CHF 1 

137 

l-chloro- 2 -hydroethforane 


s? 

Toxicitv unknown 







CCI 2 F 2 

i i 2 t 

dichlorodiiluoromethano 


s 


1 30, 51 

18,38,63 

14, 26, 63 
14, 56, 44 
14, 26 

CF 3 NF 2 

121 

perfluoromethylamine 


s? 

M usty odor. T oxicity, liquid density unknown. 

■ CCIF 3 

105 

chlorotrlfluoromethane 

Freon 13 

s 

28,30 









SOjFj-- 

102 

sulfuryl fluoride 


s? 

Toxicity, liquid density unknown. Chemistry 
resembles SFe. 

55 


14 







CHF 2 CHF 2 

102 

l-h y dro- 2 -h y droethforane 


s? 




14 


101 


Genetrou 101 

Q 

Slightly flammable, toxicity marginal. (See 
remarks for Genetron 100.) 

37 

18 




FCONF 2 

99 

carbonylnitrogentrlfluorlde 


s? 

Toxicity, liquid density unknown. Reported as 
very inert. 



14 







CF 4 -. 

88 

carbontetrafluoride 


s 

Extreme thermal stability and chemical inert- 
ness. 

28,30 

21 

14, 56, 19 





CZICIF 2 

87 

chlorodifluorometbane 


Q 

Toxicity marginal (classified U. A. L. group 5A). 

30 

38 

14 




CH 1 CF 3 

84 

1 , 1 , 1 -trifluoroethane. 


s? 


31,32 

54 


14 
14, 43 

NFj 

71 

nitrogentrifluorlde 


s? 

Toxicity uncertain since early tests (Ruff, ref. 
43) are questionable. 

43, 25 




CHFs 

70 

fluoroform. . . 


S, I 


28,30 

1G 

14 


1 



GG 

1 , 1 -difluoroc thane 

Clcnetron 100 

Q 

Slightly flammable, but forms nonflammable 
azeotropic mixture with Freon 12. Azeo- 

30 

18,25 

14 








tropic mixture called “'Carrene 7," has m= 
100, Tb= —33 ° C. 




CHjFj .... 

52 

difluoromethane.. 


s? 


30 

25 

14 







003 

44 

carbon dioxide 


Q 

Toxic for prolonged exposure at concentrations 
above about 3 per cent. 

41 

25 

40,13 





n s o 

44 

nitrous oxide 


Q 

Nonflammable but acts as oxidant at high tem- 
perature. Toxic at high concentrations. 

41 

25 

40,13 






IS satisfactory. In that sufficient Information is available to regard the gas as nontoxic, uoneorrosive, nonflammable, thermally stable to at least about 200° C, and inert chemically. 

S? may be satisfactory, but available information is insufficient to judge it so. (See remarks for each gas thus classified.) 

Q questionable because of possessing at least one undesirable characteristic (see remarks) which, however, may not be totally disqualifying. 

I inert physiotoyically, as determined by tests with small animals living in atmosphere of 20-pcreent oxygen, 80-percent gas in question, and showing no observable effects cither during 
the duration of tests (usually several hours) or afterwards. 
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TABLE III.— PHYSICAL CHARACTERISTICS OF POLYATOMIC GASES 


Formula 

7712, 

gni 

Boiling 
temp., 
T b ° C 

Critical constants 
(2) 

Viscosity 
at 0° C 
M2, 

micropoise 

?«2 

M2 

( w ) 3 /- 

Spectroscopic data fundamental wave numbers, u<, cm -1 
(3) 


mole 

(D 

7V, 

°K 

Vc, | 

atm. 

Vc 

cc/gm. mole 

CiFio 

238 

-2 

387 

23 | 

378 

111 

12.7 

0.018 


CjFdN 

221 

-11 








CA’s 

200 

-0 

390c 


33 Lc 

111 

10.7 

.023 

1431, 1385, 1340, 1285, 1239 (21, 1220, 1008 (2), 903, 745 (2), 099, GG0 (2), G13, 569, 439, 359, 










338 (2), 285, 273, 258, 250 (2), 192 (2), 173, 86 

CF,CBrF, 

199 

-22 

380c 

32c 

208c 

128 

9.2 

.031 


CFjSFs 

19(3 

— 20 








SoFs 

193 

—47 sub. 

3400 


21 Oe 

154 

7.4 

.047 

787 (3), 70S, GG2 (2), 4G1 (3), 405 (3), 245 (3) 


188 

-38 

344e, 

2Gc 

2H8e 

122 

9. 1 

.030 


C 2 F 5 X F 

171 

-37 








CBrClFs 

1(35 

-4 

427 

41.2 

232 

121 

8. 1 

.036 

1150, 1102, 872, 048, 440, 400, 380, 300, 200 

CF 3 CCIF 5 

i« 

-38 

353 

30.8 

260 

119 

7.7 

.039 

1350, 1241, 1 230, 118G, 1132, 982, 703, 648, 595, 501, 454, 441, 3GG, 331, 315, torsion, 180 (2) 

CF 3 OCF 3 

154 

-59 








CBl'F, 

149 

-60 

341 


199 

143 

G.2 

.059 

1207 (2), 1087, 702, 548 (2), 348, 297 (21 

SF» 

146 

—64 sub. 

319 

30.8 

201 

142 

6. 1 

.060 

940 (3), 775. 044 (2), 015 (3), 524 (31, 303 (3) 

CF 3 CF 3 

138 

-78 

293 

29 . n 

220 

139 

5.9 

. 063 

1420, 1250 (2), 1235 (2), 1120, 810, 715, 020 (2), 520 (2), 380 (2), 350, 214 (2), torsion 

CCIF 2 C 1 1 F 2 

137 

-10 

400 

30.7 

251 

108 

7.5 

.038 


CCl;F-.._ 

121 

-30 

385 

39.0 

218 

118 

6.1 

.055 

1159, 1101, 900, 007, 473, 455, 437, 318, 201 

CF 3 NF 2 

121 

-78 








CCIFj 

105 

-80 

302 

39 

180 

13G 

4.6 

.091 

1210 (2), 1 102, 783, 560 (2), 478, 356 (2) 

SO 2 F 2 

102 

-55 

340c 






1502, 1269, 885, 848, 553, 544, 539, 385, 300 

CHF.CHF. 

102 

-23 








CILCClFj 

101 

-9 

411 

40. 7 

232 

97 

6.2 

.049 

1 3035, 2905, 1447, 1395, 1231, 1202, 1127, 1103, 967,904, 083, 544, 526, 435, 429, 334, 305, torsion 1 

FCONKj 

91) 

-82 

1 







C Ft 

88 

-128 

228 

36. 8 

139 

171 

3. 1 

.10 

1205 (3), 904, 030 (3), 437 (2) 

CIIC1F; 

87 

-41 

370 

48, 5 

100 

120 

4.3 

.094 

1 3023, 1347, 1311, 1178, 1116, 809. 595, 422, 365 

crijCF.1 

84 

1 

1 



1 



3034 (2), 2975, 1409, 1442 (2), 1281, 1234 (2), 970 (2), 830, 606 (2), 541, 368 (2), 238, torsion 

NF, 

71 

-129 

, 2300 


125e 

104 

2. G 

.24 

1032, 647, 905 (2), 493 (2) 

CHF:i 

70 

-84 

306 

47 

136 

' 134 

3. 1 

-1G 

1370 (2), 1160 (2), 1117, 697, 508 (2), 3062 

CH 3 CHF 2 

66 

-25 

387 

44.3 

181 

90 

4.1 

.091 

3018, 3001, 2979, 2903, 1400, 1452, 1414, 1372, 1300, 1171, 1145, 1120, 925, 808, 570, 472, 392, 










torsion 

C H 2 F 2 

1 52 

I -52 

1 351 


170 

93 

3.3 

.12 

3012, 2949, 1508, 1435, 1202, 1170, 1116, 1090, 529 

co. 

44 

! —79 sub. 

1 304 

| 72.9 

96 

137 

1.9 

I .34 

2349, 1340, 067 (2) 

N>0 

44 

1 -88 

310 

71.7 

98 

135 

1 1.9 

| .33 

2237, 1288, 588 (2) 


1 sub. sublimes ul atmospheric pressure. 

2 e estimated from similar gases, or from liquid density and reduced orthobaric density of similar gases. 
3 number in parentheses denotes degree of degeneracy of wave number it follows. 


TABLE IV.— CHARACTERISTICS OF GAS-MIXTURE WIND TUNNELS 


Polyatomic gas 

Mixture propor- 
tions 

Mixture molecular weight 
relative to air 

Mixture viscosity at 0° C 
relative to air 

Size of gas-mixture wind 
tunnel relative to air wind 
tunnel for same Mm, Be, & 

Vt 

Horsepower of gas-mixture 
wind tunnel relative to air 
tunnel for same Mo,, Re, 

& pi 


7712, 

C 



argon 

krvptoii 

xenon 

argon 

krypton 

xenon 

argon 

krypton 

xenon 

argon 

krypton 

xenon 


gm. 

















mole 

mole °K 



m 

m 

m 

M 

M 

M 

L 

L 

L 

HP 

HP 

HP 



a) 





















(a) Design Mach number 1.3 

Tco =230° 

K 







CiFio 

238 

35.8c 

0.07 

0. 93 

1.82 

3.23 

4. 76 

1.13 

1.24 

1.16 

0. 84 

0.69 

0.53 

0.52 

0. 27 

0. 13 

C.Fs — 

200 

33.5 

.07 

.93 

1.77 

3. 17 

4. 69 

1.12 

1.23 

1. 15 

.84 

.09 

.53 

.53 

.27 

.13 

CF,CBtF 2 .... 

199 

24.4C 

. 10 

.90 

1.93 

3.29 

4.79 

1.12 

1,22 

1.15 

.80 

.07 

.53 

.46 

.25 

.13 

ScFo 

193 

22.5 

. 1 1 

.89 

1.98 

3.31 

4.79 

1.17 

1. 25 

1.17 

.83 

.69 

.53 

.49 

.27 

.13 

CjF s 

188 

29.5c 

.08 

.92 

1.80 

3. 18 

4.G9 

1.13 

1.24 

1. 10 

.84 

.69 

.53 

.53 

.27 

.13 

CBrClF, 

165 

15.8 

.19 

.81 

2. IS 

3.41 

4.70 

1.03 

1. 13 

1.08 

.70 

.01 

.50 

.34 

.20 

.11 

CF 3 CCIF 2 

155 

23.0 

. 11 

.89 

1.82 

3. 16 

4. 02 

1.10 

1,20 

1.13 

.81 

.08 

.53 

.49 

.20 

.13 

CBrF 3 ..__ 

149 

14.7 

.21 

.79 

2. 15 

3. 35 

4.65 

1.09 

1. 17 

1. 12 

. 74 

.64 

.52 

.37 

.22 

.13 

S Fa 

146 

18.9 

. 14 

. 80 

1.91 

3. 19 

4.59 

1. 12 

1.22 

1. 15 

.81 

.68 

.53 

.48 

.26 

.13 

CF,CF, 

138 

21.8 

.12 

.88 

1.78 

3.11 

4. 55 

1. 13 

1.23 

1.10 

.85 

.70 

.54 

.54 

.28 

.14 

CC1jF 2 

121 

15.2 

.20 

.80 

1.93 

3. 14 

4. 45 

1.02 

1. 13 

1.09 

.74 

.64 

.52 

.39 

.23 

.13 

CCIFs.. 

105 

13.8 

.23 

. / i 

1.89 

3.05 

4.35 

1.06 

1. 16 

1.12 

.77 

.67 

.54 

.44 

.26 

.14 

CH 3 CCIF 2 

101 

16.8 

. 17 

.83 

1.74 

2. 99 

4.35 

.99 

1. 12 

1.08 

.75 

.65 

.52 

.43 

.24 

.13 

CF,._ 

88 

12.5 

.27 

.73 

1.82 

2. 93 

4.14 

1.15 

1.23 

1.17 

.85 

.72 

.57 

.53 

.30 

.16 

CHCIF 2 ... .. 

87 

12.0 

.29 

.71 

1.84 

2.92 

4.10 

.98 

1.09 

1. 06 

.73 

.63 

.53 

.39 

.24 

.14 

NF, .. 

71 

11.1 

.33 

.67 

1.73 

2.74 

3. 86 

1.11 

1.20 

1. 10 

.85 

.73 

.59 

.54 

.32 

.18 

(TiF, 

70 

10.7 

.35 

.65 

1.74 

2. 72 

3. 79 

1.01 

l.n 

1.08 

.76 

.67 

.56 

.44 

.28 

.16 

CH 3 CIIF 2 

66 

13.5 

.24 

.70 

1.59 

2.74 

4.00 

.95 

1.09 

1.06 

.76 

.66 

.53 

.45 

.26 

.14 

CEIsF, 

52 

9.22 

.47 

.53 

1.57 

2.37 

3.24 

.78 

.91 

.93 

.63 

.59 

.52 

.31 

.23 

.16 

C0 2 

44 

8.11 

.63 

.37 

1.47 

2. 03 

2. 63 

.92 

1.02 

1.02 

.70 

.72 

.03 

.48 

.36 

.25 

SaO 

44 

8.45 

.57 

.43 

1.40 

2. 11 

2.81 

.94 

1.04 

1.04 

.78 

.72 

.62 

.50 

.35 

.23 







(b) Design Macli number 2.5 

Too = 140° 

K 







CBrF, 

149 

11.0 

.33 

. 67 

2. 62 

3.62 

4. 73 

1.03 

1.10 

1.06 

.64 

.58 

.49 

.25 

.17 

.11 

SFj 

146 

11.5 

.30 

.70 

2. 49 

3.54 

4.09 

1.04 

1.12 

1.07 

.66 

.59 

.50 

.28 

.19 

.11 

CF 3 CF 3 

138 

15.3 

. 19 

.81 

2. 1)3 

3. 25 

4. 57 

1.09 

1. 17 

1. 12 

.70 

.65 

.52 

.41 

.24 

.13 

CCIF 3 

105 

10.2 

.39 

.61 

2. 24 

3. 17 

4. 19 

.98 

1.07 

1.05 

.00 

.60 

.51 

.29 

.20 

.13 

CFi 

88 

9.31 

.46 

.54 

2. 14 

2.95 

3.90 

1.10 

1. 16 

1. 13 

.75 

.68 

.57 

.39 

.27 

.17 

NF, 

71 

8.74 

.53 

.47 

1.94 

2. 65 

3.43 

1.05 

1. 13 

1. 10 

.76 

.69 

.60 

.41 

.29 

.19 

CIIF, 

70 

8.60 

.55 

.45 

1.95 

2. 63 

3.38 

| .92 . 1.00 

1.00 

.00 

.02 

.54 

.31 

.23 

.16 







(c) Design Macli number 3.5 

Too = 105° 

K 







CF, 

88 

8.40 

.58 

.42 

2.34 

2.97 

3. 00 

1.07 

1.12 

1. 10 

.70 

.65 

.58 

.32 

.25 

.17 

nf 3 

71 

8.19 

.62 

.38 

2. 04 

2. 62 

3. 24 

1.03 

1.09 

1.07 

.72 

.67 

.60 

.36 

.28 

.20 


J c estimated from similar gases and from average bond frequencies of ref. 58. 
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